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Realistic nonstrange hybrid compact stars with color superconducting quark matter in their in- 
terior are constructed. It is shown that a positively charged two-flavor color superconducting phase 
could naturally appear in the core of a hybrid star as one of the components of a globally neutral 
mixed phase. The negatively charged normal quark phase is the other component of the mixed 
phase. The quark core of the star is surrounded by another mixed phase made of hadronic and 
normal quark matter. The two mixed phases are separated by a sharp interface. Finally, the lowest 
density regions of the star are made of pure hadronic matter and nuclear crust. 
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I. INTRODUCTION 

At large densities, quantum chromodynamics (QCD) 
becomes a weakly interacting theory of quarks lj. The 
quarks tend to form a highly degenerate Fermi sur- 
face. Because of the existence of an attractive interac- 
tion in the color-antitriplet channel, the famous Cooper 
instability develops, and the ground state of the sys- 
tem is rearranged 0, 0- Q • The rearrangement is sim- 
ilar to the Bardeen-Cooper-Schrieffer (BCS) mechanism 
of low temperature superconductivity in metals and al- 
loys 0. The ground state of dense quark matter is a 
color superconductor. It is characterised by spontaneous 
breaking (through the Anderson-Higgs mechanism) of 
the non-Abelian SU(3) C color gauge group rather than 
the Abelian U(l) em group of electromagnetism. 

Weakly interacting quark matter at asymptotic densi- 
ties was studied from first principles in Refs. 0, 0, @- 
These studies proved that quark matter is indeed a color 
superconductor at sufficiently high densities. An explicit 
asymptotic expression for the superconductingorder pa- 
rameter was also derived in Refs. 0, 0, IS 13- How- 
ever, the microscopic studies are not very reliable quan- 
titatively at realistic densities that exist in nature (e.g., 
inside compact stars). The corresponding running cou- 
pling constant of QCD defined at the relevant scale of 
the quark chemical potential is large, and therefore the 
results cannot be trusted. 

It might be appropriate to mention that the creation 
of cold dense quark matter with color superconducting 
properties in heavy ion collisions seems very unlikely to 
us (see, however, some speculations in Ref. |lOj). In par- 
ticular, it does not look plausible that one could avoid 
producing a large entropy per baryon in any type of a 
heavy ion collision. Therefore, a cold and dense environ- 
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ment, needed to support color superconductivity, could 
hardly be formed. A typical central region of a compact 
star, on the other hand, appears to be a very natural 
place where sufficiently cold and dense matter could ex- 
ist. Indeed, the estimated central densities of such stars 
might be as large as 10po (where po » 0.15 fm -3 is the 
saturation density), while their temperatures are in the 
range of 10 to 100 keV. 

In order to study quark matter at realistic densities, 
currently there is no other real alternative but to rely on 
various phenomenological models (such, for example, as 
Nambu- Jona-Lasinio type models 0, El O El ) ■ The 
obvious shortcoming of this approach is the same as that 
of the microscopic approach: the results should always 
be interpreted with caution and treated as qualitative 
rather than quantitative. Nevertheless, if one accepts the 
possibility of color superconductivity in quark matter at 
realistic densities, existing in compact stars, it is not so 
important which specific model is used to estimate the 
effect of color superconductivity on properties of compact 
stars. 

In general, QCD at high baryon density has a very 
rich phase structure. There are many possible color su- 
perconducting phases of quark matter made of one, two 
and three lightest quark flavors. Each of them is charac- 
terised by a unique symmetry breaking pattern and by 
a specific number of bosonic as well as fermionic gapless 
modes. In the rest of this paper, we are going to con- 
centrate our attention almost exclusively on matter with 
two quark flavors. In order to better understand our mo- 
tivation, let us also briefly mention other possibilities. 

A one-flavor color superconductor is characterised by 
a condensate of Cooper pairs made of the same quark 
flavor. Thus, the corresponding wave function is sym- 
metric in flavor. By taking into account that the con- 
densation results from the attractive interaction between 
quarks in the color antitriplet antisymmetric channel, one 
finds that a spin zero configuration is forbidden by the 
Pauli principle. Therefore, it is a much weaker spin-1 
condensate that gives rise to color superconductivity in 
one-flavor dense quark matter The recent study 

in Ref. 0] suggests that the presence of such a phase 
may still affect some properties of compact stars. 
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Perhaps, the most interesting color superconducting 
phase appears in quark matter with three light flavors. 
It is the so called color-flavor-locked (CFL) phase (l2^ . 
In this phase the original gauge symmetry SU(3) C and 
the global chiral symmetry SU(3) l X SU(3)# break down 
to a global diagonal SU(3) c +l+r subgroup. As a re- 
sult, the quark quasiparticles get large gaps in their en- 
ergy spectra and decouple from the low energy dynam- 
ics. Because of the Higgs mechanism all gluons also be- 
come massive and decouple. Then, it appears that the 
low energy spectrum of the CFL phase is dominated by 
one Nambu-Goldstone boson and nine massive pseudo- 
Nambu-Goldstone bosons This observation has in- 
teresting implications for the transport properties of the 
CFL phase 0. 

In order to have the CFL phase in compact stars, the 
strange quark should be sufficiently light at the corre- 
sponding densities. However, the actual density depen- 
dence of the constituent, medium modified mass of the 
strange quark is not known in QCD. Taking this into ac- 
count, one could imagine that the strange quark might 
be too heavy even to appear in baryonic matter at realis- 
tic conditions inside compact stars. Then, of course, the 
CFL quark phase could not be realized. In yet another 
scenario, the strange quarks could appear inside stars, 
but they would not participate in Cooper pairing. Then, 
a mixture of 2SC matter and normal matter of strange 
quarks could form. We denote such a phase 2SC+S. 

Matter in the bulk of a compact star should be neutral 
(at least, on average) with respect to electric as well as 
color charges. Otherwise, the gravitational force would 
not be able to hold the star together. Also such mat- 
ter should remain in /3-equilibrium. These conditions 
turn out to play an important role in determining which 
color superconducting phases can and which cannot exist 
inside stars. Indeed, the strongest color superconduct- 
ing phases appear when the Fermi momenta of differ- 
ent quark flavors are approximately equal. By enforcing 
the neutrality and /3-equilibrium conditions, however, one 
does not necessarily get the Fermi surfaces of quarks that 
are best suitable for Cooper pairing. In some cases, these 
conditions may even appear to be inconsistent with color 
superconductivity. 

Recently, it was argued in Ref. 01 that imposing the 
charge neutrality condition on quark matter inside com- 
pact stars would prevent the appearance of the 2SC+S 
phase, and would favor the appearance of the CFL phase 
instead. (Note that the strange quark mass was chosen 
too small to allow the appearance of the pure 2SC phase 
in Ref. H3-) This conclusion was essentially confirmed 
m a more rigorous study of Ref. [3 where it was found 
that the 2SC+S phase could exist only in a narrow win- 
dow (about 10 or 15 MeV wide) of the baryon chemical 
potential around the midpoint /ib/3~ 450 MeV. In this 
phase, a positive charge of the 2SC condensate is compen- 
sated by a negative charge of unpaired strange quarks. At 
smaller values of the chemical potential, no conventional 
2SC was detected in Ref. 0|. It was shown in Ref. 0] 



that a new type of neutral 2SC phase, the so called gap- 
less 2SC phase develops at the corresponding densities. 
In this gapless 2SC phase, the conventional relations be- 
tween the number densities of the pairing quarks are not 
valid. 

While matter in the bulk of a star should satisfy the 
condition of charge neutrality, this condition should not 
necessarily be enforced locally. In fact, it is sufficient to 
make matter neutral on average, or globally. This kind 
of matter appears naturally in mixed phases. Recently, 
studies of various mixed phases with color superconduct- 
ing matter were presented in Ref. pi). In this paper, 
we use this general idea to construct a realistic model 
of a compact star with two-flavor color superconducting 
matter at the core (for constructions using_2SC quark 
matter without mixed phases see Ref. [22|, |23| , and for 
constructions using CFL quark matter see Ref. pi . [2j| ) . 

This paper is organized as follows. In the next section, 
we present a quark model with a four-fermion interac- 
tion that we use to derive the quark matter equation of 
state. There, we also briefly discuss the phase diagram of 
the model, and the subtleties of enforcing the electrical 
neutrality condition in color superconducting matter. In 
Sec. IIIII we present the Gibbs construction for the mixed 
phase of normal and color-superconducting quark matter. 
In Sec. IIVI we discuss a model for hadronic matter, and 
the corresponding equation of state with the charge neu- 
trality condition imposed. In the same section, we apply 
the Gibbs construction to derive the equation of state of 
the hadron-quark mixed phase. The star properties are 
discussed in Sec. and the summary and outlook are 
given in Sec. IVII 

II. QUARK MODEL 

We start our analysis from discussing the quark model 
that we use to derive the equation of state of baryonic 
matter at high density. In this study, we accept a rather 
conservative point of view that the strange quark is suf- 
ficiently heavy and does not appear at baryon densities 
that can be reached in compact stars without causing a 
gravitational instability. It is also appropriate to men- 
tion that we aim at a general qualitative rather than 
quantitative description of hybrid stars with color super- 
conducting matter in their interior. These assumptions 
justify the use of the simplest SU(2) Nambu-Jona-Lasinio 
model of Ref. p(j in our study. The explicit form of the 
Lagrangian density reads: 

C = q(i'fdp - m )q + G s [{qq) 2 + [qi^Tqf] 

+ G D [(iq C ee h l5 q){zqe^l,q C )\ , (1) 

where q c — Cq T is the charge-conjugate spinor and 
C = ij 2 j° is the charge conjugation matrix. The quark 
field q = qi a is a four-component Dirac spinor that car- 
ries flavor (i = 1,2) and color (a = 1,2,3) indices. 
t = (t 1 , t 2 , t 3 ) are Pauli matrices in the flavor space, 



3 



and (e) 



; k 



Ma? 



f af3b 



are antisymmetric ten- 



sors in flavor and color, respectively. We also introduced 
two independent coupling constants in the scalar quark- 
antiquark and scalar diquark channels, Gs and Gd- The 
definition of this non-renormalizable NJL model is com- 
plete after introducing a momentum cut-off A. 

The values of the par ameters in the NJL model are 
the same as in Refs. [HIU: G s = 5.0163 GeV" 2 and 
A = 0.6533 GeV. Here we also consider only the chiral 
limit with m = 0. As for the strength of the diquark 
coupling Gd, its value is taken to be proportional to 
the quark-antiquark coupling constant, i.e., Gd = V^S 
with rj = 0.75. The choice 77 = 0.75 corresponds to the 
regime of intermediate strength of the diquark coupling 
|l9| . Regarding this point, we mention that our final 
results would remain qualitatively the same also for a 
wide range of weak and strong diquark couplings. 

In /3-equilibrium, the diagonal matrix of quark chemi- 
cal potentials is given in terms of baryonic, electric and 
color chemical potentials. In particular, 



One should notice that, in general, there are two genera- 
tors (T 3 and T 8 ) in the center of the SU(3) C color group. 
Therefore, one could introduce two chemical potentials 
for two different color charges. However, we require that 
quark matter in the 2SC ground state remain invariant 
under the SU(2) C color gauge subgroup. This condition 
forbids the introduction of the second nontrivial color 
chemical potential /J3. 

The effective potential for quark matter at zero tem- 
perature and in /3-equilibrium with electrons takes the 
form [H 
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where fto is a constant added to make the pressure of the 
vacuum zero. The sum in Eq. runs over all (6 quark 
and 6 antiquark) quasiparticles. The dispersion relations 
and the degeneracy factors of the quasiparticles read 



Hij,a0 = (fJ-Sij — ^eQij)S a /3 H -r= flS^ij (Tg) a /3 , (2) 

v3 

where Q and T 8 are generators of U(l) em of electromag- 
netism and U(l)s subgroup of the color gauge group. 
The explicit expressions for the quark chemical poten- 
tials read 
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^/[E{p)±p} 2 + A 2 ±6^. [x2] (10) 
Here we used the following shorthand notation: 



E* = E(p)±fM ub , 
E{p) ± fi db , 
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After making use of the explicit quasiparticle dispersion 
relations, we arrive at the following potential: 
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where p^ = p ± yj (Sfi) 2 — A 2 . Note that the physical 
thermodynamic potential that determines the pressure, 
ftphys = —P, is obtained from ft in Eq. l|14f) after sub- 
stituting the order parameter A that solves the following 
gap equation: 



9ft „ 



(15) 



In addition, homogeneous quark matter in the bulk of a 
compact star should be neutral with respect to color and 



electric charges. Therefore, for such matter, one should 
also impose extra two (local) conditions: 
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The solution that satisfies the gap equation and both 
neutrality conditions was studied in detail in Ref. [l9| . It 
was found that the corresponding phase of matter is the 
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gapless 2SC phase. This phase has the same symmetry 
of the ground state as the conventional 2SC phase. In 
the low-energy spectrum, however, it has two additional 
gapless fermionic quasiparticles. 

In order to determine the most favorable phase of 
quark matter inside a hybrid compact star, one has to 
consider all possible neutral phases, and compare their 
values of the pressure at a given baryon chemical poten- 
tial. The simplest candidates are the gapless 2SC and 
normal quark phases. Both of these can be made locally 
neutral. For any given value of the baryon chemical po- 
tential, gapless 2SC matter has a slightly higher pressure 
than normal quark matter |l9j . Numerically, we find 
that the pressure difference SP lies somewhere between 
1 and 2.5 MeV/fm 3 . Therefore, the gapless 2SC phase 
is slightly more favorable than the neutral normal quark 
phase. 

In addition to homogeneous (one-component) phases, 
one could also study various mixed phases of two-flavor 
quark matter. The most promising components for con- 
structing mixed phases are (i) normal quark matter, (ii) 
gapless 2SC matter, and (iii) ordinary 2SC matter. The 
first two of them are included because they allow locally 
neutral phases by themselves. The last one, on the other 
hand, is included because it is the strongest color su- 
perconducting phase known. Clearly, this list of phases 
is not complete. One could also take into consideration 
other phases, e.g., various combinations of color super- 
conductors with spin-1 condensates 0,EIE|- Our expe- 
rience suggests that these latter would lead to less favor- 
able one- and multi-component constructions. In order to 
prove this rigorously, a detailed study is required. This, 
however, is outside the scope of the present paper. 

Inside mixed phases, the charge neutrality is satisfied 
"on average" rather than locally. This means that dif- 
ferent components of mixed phases may have non-zero 
densities of conserved charges, but the total charge of all 
components still vanishes. In this case, one says that the 
local charge neutrality condition is replaced by a global 
one. Now, in absence of the local neutrality, the pres- 
sure of each phase could be considered as a function of 
baryon chemical potential, as well as a function of chemi- 
cal potentials related to other conserved charges (e.g., /x e 
and fj,g in the case of two- flavor quark matter at hand). 
In order to decide which mixed phase is the most favor- 
able, one should utilize the so-called Gibbs construction. 
The details of this construction will be given in the next 
section. 

While we intend to relax the local neutrality condition 
with respect to the electric charge, in this paper we al- 
ways enforce the condition of local color neutrality. We 
cannot fully justify this constrain in all situations where 
mixed phases appear. However, one could speculate that 
color charged components are less likely to make ener- 
getically favorable mixed phases because of strong color 
Coulomb forces. Besides, in the two main constructions 
that we build in this paper, at least one of the com- 
ponents (normal quark or hadronic matter) of a mixed 




FIG. 1: Pressure as a function of /j, = and /i e for the 

normal and color superconducting quark phases. The dark 
solid lines represent two locally neutral phases: (i) the neutral 
normal quark phase on the left, and (ii) the neutral gapless 
2SC phase on the right. The appearance of the swallowtail 
structure is related to the first order type of the phase tran- 
sition in quark matter. 



phase is expected to be locally color neutral. Therefore, 
our consideration here remains rather general. 

The pressure of the main three phases of two-flavor 
quark matter as a function of the baryon and electric 
chemical potentials is shown in Fig. ^ In this figure, 
we also show the pressure of the neutral normal quark 
and gapless 2SC phases (two dark solid lines). In full 
agreement with the study of Ref. [l9| , the surface of the 
gapless 2SC phase extends only over a finite range of the 
values of /i e . It merges with the pressure surfaces of the 
normal quark phase (on the left) and with the ordinary 
2SC phase (on the right). 

It is interesting to notice that the three pressure sur- 
faces in Fig. n form a characteristic swallowtail struc- 
ture. As one could see, the appearance of this structure 
is directly related to the fact that the phase transition 
between color superconducting and normal quark mat- 
ter, which is driven by changing parameter fi e , is of first 
order. In fact, one should expect the appearance of a 
similar swallowtail structure also in a self-consistent de- 
scription of the hadron-quark phase transition. Such a 
description, however, is not available yet. 

^From Fig. ^ one could see that the surfaces of normal 
and 2SC quark matter intersect along a common line. 
This means that the two phases have the same pressure 
along this line, and therefore could potentially co-exist. 
Moreover, as is easy to check, normal quark matter is 
negatively charged, while 2SC quark matter is positively 
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charged on this line. This observation suggests that the 
appearance of the corresponding mixed phase is almost 
inevitable. The details of the constructions are given in 
the next section. 




III. GIBBS CONSTRUCTION WITHIN QUARK 
MODEL 



In this section we study the possibility of mixed phases 
in two-flavor quark matter. As we have already men- 
tioned in the previous section, we restrict our consid- 
eration to only three most promising components: (i) 
normal quark matter, (ii) gapless 2SC matter, and (hi) 
ordinary 2SC matter. Before proceeding further, it would 
be appropriate to mention that some mixed phase with 
color superconducting components were also constructed 
in Ref. |2l|. 

Let us start by giving a brief introduction into the 
general method of constructing mixed phases by impos- 
ing the Gibbs conditions of equilibrium [2H 1271 ] . From 
the physical point of view, the Gibbs conditions enforce 
the mechanical as well as chemical equilibrium between 
different components of a mixed phase. This is achieved 
by requiring that the pressure of different components 
inside the mixed phase are equal, and that the chemi- 
cal potentials (fj, and ji e ) are the same across the whole 
mixed phase. For example, in relation to the mixed phase 
of normal and 2SC quark matter, these conditions read 



/' 

/'<: 



p(2SC), 
(NQ) 



(2SC) 
(2SC) 



(18) 
(19) 
(20) 



It is easy to visualize these conditions by plotting the 
pressure as a function of chemical potentials (/z and /z e ) 
for both components of the mixed phase. This is shown 
in Fig. [21 As should be clear, the above Gibbs conditions 
are automatically satisfied along the intersection line of 
two pressure surfaces (dark solid line in Fig. |2J). 

Different components of the mixed phase occupy differ- 
ent volumes of space. To describe this quantitatively, we 
introduce the volume fraction of normal quark matter 

as follows: X2.sc = Vnq/V (notation \b means "vol- 
ume fraction of phase A in a mixture with phase B"). 
Then, the volume fraction of the 2SC phase is given by 
xWq — (1 — X2Sc)- F rom t ne definition, it is clear that 

0<X^c<l- ^ 

The average electric charge density of the mixed phase 
is determined by the charge densities of its components 
taken in the proportion of the corresponding volume frac- 
tions. Thus, 

n[ MP) = X^C^VMe) + (1 - x^>i 2SC) (^e). 

( 21 ) 

If the charge densities of the two components have op- 
posite signs, one can impose the global charge neutrality 



0.25 




0.2 



0.15 



0.1 



0.05 



FIG. 2: Pressure as a function of n = and fi e for the 

normal and color superconducting quark phases (the same as 
in Fig. but from a different viewpoint). The dark solid 
line represents the mixed phase of negatively charged normal 
quark matter and positively charged 2SC matter. 



could not exist. In the case of quark matter, the charge 
density of the normal quark phase is negative, while the 
charge density of the 2SC phase is positive along the line 
of the Gibbs construction (dark solid line in Fig. 
Therefore, a neutral mixed phase exists. The volume 
fractions of its components are 



NQ 
X2SC 



V 2SC 

Xnq 
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(2SC) (NQ) ' 
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, NQ n 

1 - X 2S C = 



(NQ) 
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(22) 
(23) 



After the volume fractions have been determined from 
the condition of the global charge neutrality, we could 
also calculate the energy density of the corresponding 
mixed phase, 



s {Mp) =x N 2 S y [NQ) ^^) 



(1 



condition, n 



(MP) 



0. Otherwise, a neutral mixed phase 



This is essentially all that we need in order to construct 
the equation of state of the mixed phase. 

So far, we were neglecting the effects of the Coulomb 
forces and the surface tension between different compo- 
nents of the mixed phase. In a real system, however, 
these might be important. In particular, the balance be- 
tween the Coulomb forces and the surface tension deter- 
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mines the geometries of different components inside the 
mixed phase. 

In our case, nearly equal volume fractions of the two 
quark phases are likely to form alternating layers (slabs) 
of matter. The energy cost per unit volume to produce 

such layers scales as cr 2 / 3 (ni 2SC ^ — ni^-*) 2 / 3 where a 
is the surface tension [2^]. Therefore, the quark mixed 
phase is a favorable phase of matter only if the surface 
tension is not too large. Our simple estimates show that 
o'max ^ 20 MeV/fm 2 . However, even for slightly larger 
values, 20 < a < 50 MeV/fm 2 , the mixed phase is still 
possible, but its first appearance would occur at larger 
densities, 3po < Pb 5po- The value of the maximum 
surface tension obtained here is comparable to the esti- 
mate in the case of the hadronic-CFL mixed phase ob- 
tained in Ref. |2||. The thickness of the layers scales as 

( T 1 / 3 (ni 2SC) - ni NQ) )- 2 / 3 H, and its typical value is of 
order 10 fm in the quark mixed phase. This is similar 
to the estimates in various hadron-quark and hadron- 
hadron mixed phases j2^, |2j| . While the actual value of 
the surface tension in quark matter is not known, in this 
study we assume that it is not very large. Otherwise, 
the homogeneous gapless 2SC phase should be the most 
favorable phase of nonstrange quark matter |19| . 

Under the assumptions of this paper, the mixed phase 
of normal and 2SC quark matter is the most favorable 
neutral phase of matter in the model at hand. This 
should be clear from observing the pressure surfaces in 
Figs. n an d0 For a given value of the baryon chemical 
potential p — the mixed phase is more favorable 

than the gapless 2SC phase, while the gapless 2SC phase 
is more favorable than the neutral normal quark phase. 

The validity of the quark model is limited when the 
baryon density decreases. From the physical point of 
view, at some point quark matter should become con- 
fined and the NJL model should fail to reproduce the 
correct equation of state of low density baryonic mat- 
ter. Therefore, at low densities, it is appropriate to use a 
hadronic description of matter. This is discussed in the 
next section. 



IV. LOW DENSITY REGION AND 
QUARK-HADRON PHASE TRANSITION 

At densities around normal nuclear matter density 
Po ~ 0.15 fm -3 , the description of baryonic matter in 
terms of quarks could hardly be adequate. At such low 
densities, quarks are confined inside hadrons. Thus, it is 
more natural to use an effective hadronic model. 

In the literature, there exist many hadronic models 
that have been studied in detail. In particular, the 
Walecka model [3(3, |^| and its nonlinear extensions 
[32I l33l | are well known. These models have been quite 
successful and widely used for the description of finite nu- 
clei, hadronic and neutron star matter. One could also 
use some hadronic versions of the NJL model [341 l35l l36| . 
as well as many others [13, HU . 



In the following we use a QCD-motivated hadronic 
chiral SU(3)lxSU(3),r model as an effective theory of 
strong interactions to describe the low density regime of 
the baryonic matter j3^,^3>^|. This model was found to 
describe reasonably well the hadronic masses of various 
SU(3) multiplets, finite nuclei, hypernuclei, excited nu- 
clear matter and neutron star properties [3^.l4fll4l| . The 
basic features of the chiral model are: (i) Lagrangian of 
the model is constructed in accordance with the nonlinear 
realization of the chiral SU(3)lxSU(3)h symmetry; (ii) 
heavy baryons and mesons get their masses as a result 
of spontaneous symmetry breaking; (iii) the masses of 
pseudoscalar mesons (pseudo-Nambu-Goldstone bosons) 
result from an explicit symmetry breaking; (iv) a QCD- 
motivated field that describes the gluon condensate (dila- 
ton field) is introduced. 

It is expected that the phase transition between the 
hadronic phase and the normal quark phase is a first 
order phase transition at zero temperature and finite 
baryon chemical potential. Then, the hadronic and quark 
phases could co-exist in a mixed phase [2^, |2]|. This 
mixed phase should satisfy the Gibbs conditions of equi- 
librium which are similar to those discussed in the previ- 
ous section, see Eqs. (tTHll - <{2"0")) . The total energy density 
in the hadron-quark mixed phase is given by 

£ (mp) = x £Q e 0vQ) (M)/ie) + (1 _ x £Q )e (*) (M)/ie); (25) 

where Xh^ denotes the volume fraction of normal quark 
matter inside a mixture with hadronic matter. 

To visualize the Gibbs construction of the hadron- 
quark mixed phase, we plot the hadronic surface of the 
pressure along with the quark surfaces discussed in the 
previous section. Thus, Fig. |2Jis replaced by Fig. [21 The 
new figure shows the surfaces of the pressure of the pure 
hadronic and quark phases as a function of their chemi- 
cal potentials p and p e . The intersection lines of differ- 
ent surfaces indicate all potentially viable mixed phase 
constructions. Although the Gibbs conditions are satis- 
fied along all these lines, not all of them could produce 
globally neutral phases (e.g., there are no neutral con- 
structions along the light shaded solid lines in Fig.[3J). 

The dark solid line gives the complete, most favor- 
able construction of globally neutral matter in general 
/3-equilibrium. This line consists of three pieces. The 
lowest part of the curve (up to the point denoted by □ 
and P < 10 MeV/fm 3 ) corresponds to the pure confined 
hadronic phase. Within this region matter is mostly com- 
posed of neutrons with little fractions of protons and elec- 
trons to realize the charge neutrality and /3-equilibrium. 
The hyperonic particles (A, £ and S) are not present in 
this lowest density region. Such particles would appear 
in the hadronic phase at considerably higher values of 
pressure and density. 

The mixed phase of hadronic and normal quark mat- 
ter starts at the baryonic density ps ~ l-49po which 
corresponds to the D-point in Fig. [31 At this point the 
first bubbles of deconfined quark matter appear in the 
system. At the beginning of this hadron-quark mixed 
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FIG. 3: Pressure as a function of fi = /j.b/3 and fi e for the 
hadronic phase (at the bottom), for the two- flavor color su- 
perconducting phase (on the right in front) and the normal 
phase of quark matter (on the left, and back on the right). 
The dark thick line follows the neutrality line in hadronic mat- 
ter, and two mixed phases: (i) the mixed phase of hadronic 
and normal quark matter; and (ii) the mixed phase of normal 
and color superconducting quark matter. 

phase, the deconfined bubbles are small but highly neg- 
atively charged, whereas hadronic matter, in which the 
bubbles are embedded, is slightly positively charged. The 
global charge neutrality condition reads 

EE X^^We) + (1 - X £Vi H We) = 0. 

(26) 

where ni H ' and Tie are the charge densities of hadro- 
nic and normal quark matter, respectively. This condi- 
tion should be satisfied at each point along the middle 
part of the dark solid line (i.e., between the points de- 
noted by □ and A). With increasing density (from about 
1.49;Oo up to 2.56po)j the volume fraction of hadronic 
matter decreases (down to about 0.59), while the fraction 
of normal quark matter increases (up to about 0.41). 

It is quite remarkable that the hadron-quark mixed 
phase does not reach a point where the fraction of 
hadronic matter would vanish completely. This is in con- 
trast to other examples of hadron-quark Gibbs construc- 
tions in the literature [2rl l42l |43| . Instead, at baryon 
density about 2.56po (a point denoted by A in Fig.[3Jl, the 
mixed phase is replaced by another mixed phase which 
is made of the normal and 2SC quark components. At 
this point, positively charged hadronic matter will be 
suddenly converted into positively charged color super- 
conducting quark matter in the 2SC phase. As a re- 
sult of this rearrangement, the values of the baryon den- 
sity and the energy density experience small jumps: the 
baryon density changes from about 2.56po to 2.75po and 




e[GeV/fm 3 ] 

FIG. 4: The equation of state for globally neutral hybrid 
matter (solid line) and globally neutral quark matter (dashed 
line). The points of the beginning of the two mixed phases 
are denoted by □ and A. 



the energy density changes from 378 MeV/fm 3 to 415 
MeV/fm 3 . Right after the transition, the volume frac- 
tions of the 2SC and normal quark phases are 0.53 and 
0.47, respectively. 

The mixed phase of normal and 2SC quark matter 
above the point A in Fig. [21 will remain the most favorable 
globally neutral phase of baryonic matter in our model. 
The volume fractions of the two components of this phase 
stay nearly constant with increasing density. This un- 
usual property of the mixed phase originate from the sim- 
ple fact that 2SC quark matter has a strong preference 
to remain positively charged. In principle, such a mixed 
phase could possibly be replaced by a neutral 2SC+S or 
CFL phase (or, possibly, even by some new mixed phase) 
at sufficiently high densities when the strange quarks ap- 
pear in the system. However, it may also happen that 
the stars with high enough densities become gravitation- 
ally unstable long before the strange quarks get a chance 
to appear. 

The equations of state for quark and hybrid matter are 
shown in Fig. The first equation of state corresponds 
to globally neutral quark matter which is a mixture of 
the normal quark and 2SC phases. This was discussed 
in Sec. IIIII in detail, see Fig. [21 As for the equation of 
state of hybrid matter, it is constructed out of the equa- 
tion of state of neutral hadronic matter and two Gibbs 
constructions in accordance with Fig. [3J As before, the 
points that indicate the beginning of two different mixed 
phases are denoted by □ and A in Fig. 
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V. STAR PROPERTIES 

It is evident from Fig. 0] that no minimum appears 
in the equation of state of hybrid matter. Therefore, 
the corresponding compact stars can only be bound by 
the gravitational force. Because of a strong gravitational 
field inside and around such compact stars, we need to 
describe the gravitational field within the framework of 
general relativity as a curvature of spacetime. By know- 
ing the behavior of the pressure and the energy density of 
matter inside the star, one can obtain the corresponding 
spacetime geometry by solving the Einstein equation, 



(27) 



where is the Ricci tensor, g^ v is the spacetime metric 
and Tn U is the energy-momentum tensor of matter inside 
the compact star. 

The metric inside a non-rotating spherically symmet- 
ric star with the energy-momentum tensor T^„ of an 
ideal fluid |44| can be determined by solving the Tolman- 
Oppenheimer-Volkoff (TOV) equations [4{| that follow 
from the above Einstein equation. The metric is given 
by the following ansatz: 



— r 2 , — r 2 sin 2 £ 



.^=diag(W>,-(l-^y\ 

and the explicit form of the TOV equations read 

m(r) = 4tt j e{r\)r\dri 1 

dv 
dr 
dP 
dr 



m(r) + Airr 3 P(r) 
r [r — 2m(r)] 



= -[P(r) + e(r)] 



dv 
dr 



(28) 

(29) 
(30) 
(31) 



(To simplify notation, we use units with k = c = 1.) 
The functions P(r) and e(r) are the pressure and energy 
density at radius r inside the star. By definition, R is the 
radius of the star while M = m(R) is the star mass. For 
a given equation of state P(e) and a fixed central energy 
density e c = e(r = 0) one can numerically integrate the 
above set of equations from the centre of the star up to its 
surface (r = R) where the pressure is zero, i.e., P(R) = 0, 
and 



. . 1 / 2M 
v{r)= -In II - — 



for r > R. 



(32) 



The energy density profiles for hybrid stars with different 
central energy densities are displayed in Fig. [S] The star 
with the lowest central energy density in Fig. [21 e c = 210 
MeV/fm 3 , is composed of pure confined hadronic matter, 
mainly neutrons, surrounded by a thin compact star crust 
consisting of leptons and nuclei. To get the equation of 
state for the crust, we use the results of Ref. |4o| for 
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FIG. 5: Energy density profiles for hybrid stars. The loca- 
tions of the interface between the two types of mixed phases 
are denoted by A, while the locations of the boundary be- 
tween the pure hadronic phase and the hadron-quark mixed 
phase are denoted by □. 



p B < 0.001 far 3 and the results of Ref. gg for 0.001 
far 3 < p B < 0.08 far 3 . 

The next energy density profile in Fig. [5] corresponds 
to the central energy density e c = 370 MeV/fm 3 . We see 
that the corresponding star already has a rather large 
core (the radius is about 8 km) consisting of a mixture of 
hadronic and normal quark matter. The core of the star 
is surrounded by a layer of hadronic matter and a crust. 

As we saw in the previous section, there is no globally 
neutral baryonic matter that would produce the energy 
densities in the window between 378 MeV/fm 3 and 415 
MeV/fm 3 . Therefore, there are no stars that could have 
the central energy densities in this window either. At 
e c > 415 MeV/fm 3 , a quark core (made of the mixed 
phase of normal quark and 2SC matter) forms at the 
center of the star. Two examples of the corresponding 
energy density profiles are also shown in Fig. |SJ The 
star with the central energy density e c = 500 MeV/fm 3 
contains a quark phase core with radius about 6 km. 
This core is separated from the layer of the hadron- 
quark mixed phase by a sharp interface (the correspond- 
ing point is denoted by A in Fig. [SJ. The most extreme 
stable star within this model (the star with the largest 
possible mass), has the central energy density e c = 1392 
MeV/fm 3 . This star is mainly composed of the quark 
core surrounded by a relatively thin layer of the hadron- 
quark mixed phase, as well as the pure hadronic phase 
and the crust on the outside. 

The dependence of the gravitational mass M as a func- 
tion of the central energy density e c for hybrid and quark 
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FIG. 6: The gravitational mass M versus the central energy 
density for hybrid stars (solid line) and quark stars (dashed 
line). The corresponding equations of state are shown in Fig. 
|1| The stars heavier than the "O-star" have central baryon 
densities larger than 5po- Such stars may have seeds of strange 
quark matter at their cores. 



stars are displayed in Fig. EI The first appearance of the 
two mixed phases inside hybrid stars are marked with 
the symbols □ and A. The maximum mass hybrid star 
has the following properties: M max = 1.81M©, e c = 1392 
MeV/fm 3 , p c = 7.58po and R = 10.86 km. Hybrid stars 
with higher masses or central densities would collapse to 
black holes because of gravitational instability. 

The mass-radius relations for hybrid and quark stars 
are shown in Fig. [7] As one would expect, the pure quark 
stars composed of the mixed phase constructed in Sec. IIIII 
have much smaller radii and the value of their maximum 
mass is slightly smaller, see Fig. [7| The difference be- 
tween hybrid and pure quark stars is mostly due to the 
low density part of the equation of state. This is also 
evident from the qualitative difference in the dependence 
of the radius as a function of mass for the hybrid and 
quark stars with low masses. The corresponding hybrid 
stars are large because of a sizable low density hadronic 
layer, while the quark stars are small because that have 
no such layers. 

Our results for pure quark stars are comparable to 
those in Refs. jH |H IH H^. Also, the maximum 
masses and the corresponding radii of the hybrid stars 
obtained here are similar to those of the strange hybrid 
stars of Ref. > assuming that the strange quark mass 
is not very small (m s > 300 MeV) and the supercon- 
ducting gap is not too large (A < 50 MeV). At smaller 
values of the strange quark mass and/or larger values of 
the superconducting gap, the strange hybrid stars tend 



FIG. 7: The mass-radius relations for hybrid stars (solid line) 
and quark stars (dashed line). 

to have smaller maximum masses and smaller radii |24j . 

In this paper we use a two-flavor version of the NJL 
model to describe the deconfined phase. However, the 
three-flavor extension of the NJL model in Ref. 0] 
suggests that strange quarks might be present in mat- 
ter above a critical density of about 5po. In accordance 
with this, we have marked the points of the expected 
appearance of strange quarks in the center of the stars 
with the symbol O i n Figs. H3 and [7| For a more realistic 
picture of hybrid and quark stars with masses and cen- 
tral energy densities larger than those of the "O-star" , 
strange quarks may need to be included in the descrip- 
tion. The corresponding study of compact stars is left for 
future work. It is worth to emphasize, however, that the 
properties of the hybrid stars lighter than the "O-star" 
in Figs. and would remain unchanged. 

VI. SUMMARY AND OUTLOOK 

In this paper, we constructed a realistic equation of 
state of nonstrange baryonic matter that is globally neu- 
tral and satisfies the condition of /3-equilibrium. This 
equation of state might be valid up to densities of about 
bpo, or even up to densities of about 8po in the most 
optimistic scenario. In our construction, matter at low 
density (ps < 1.49po) is mostly made of neutrons with 
traces of protons and electrons. At intermediate densi- 
ties (1.49po ^ Pb ^ 2.56po), homogeneous hadronic mat- 
ter is replaced by the mixed phase of positively charged 
hadronic matter and deconfined bubbles of negatively 
charged normal quark matter. The volume fraction of 
normal quark matter gradually grows with increasing 
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density. Before the volume fractions of two phases be- 
come equal, the mixed phase undergoes a rearrangement 
in which the hadronic component of matter turns into 
the two-flavor color superconductor. At higher densities 
(ps > 2.75po)> only the quark mixed phase exists. This 
latter is composed of about equal fractions of normal and 
2SC quark matter. 

Previously it was argued that 2SC quark matter could 
not appear in compact stars when the charge neutrality 
condition is imposed |17| . The main reason for this is 
a strong preference of 2SC quark matter to remain pos- 
itively charged. Therefore, one may conclude that no 
color superconducting phases can appear inside compact 
stars if the baryon density is not large enough for the 
strange quarks to condense and participate in the for- 
mation of the neutral CFL phase. The present study 
reinstalls the status of 2SC quark matter as the most 
promising color superconducting phase in the central re- 
gions of compact stars. 

Our study shows that positively charged 2SC quark 
matter appears naturally in a quark mixed phase at den- 
sities around 3po- The other component of the globally 
neutral mixed phase is negatively charged normal quark 
matter. The corresponding construction turns out to be 
very stable. In particular, we observe that the volume 
fractions of the two quark components remain approxi- 
mately the same with changing the baryon density. From 
the physical point of view, this rigidity of the mixed phase 
is connected with the nature of 2SC quark matter which 
tends to remain positively charged. 

To the best of our knowledge, hybrid baryonic mat- 
ter in this paper gives the first and the only example of 
a combination of two different Gibbs constructions that 
replace each other in the same system with changing the 
density. It should be clear that this special construc- 
tion results from the existence of the triple point where 
hadronic matter coexists with two different quark phases. 
It is fare to mention that some speculations about the 
possibility of another type of a triple point were expressed 
in Ref. j4£|. In our construction, at the triple point, the 
hadronic component of the first mixed phase is replaced 
by the 2SC quark component of the other. This could 
be interpreted as a first order phase transition that hap- 
pens only in one of the components of the mixed phase. 
The volume fraction of the inert component, i.e., nor- 
mal quark matter, changes by a small jump during this 
transition. This change accounts for the difference of the 
charge densities of hadronic and 2SC quark matter. 

By making use of the equation of state of hybrid mat- 
ter, we construct the corresponding non-rotating com- 
pact stars. We find that the largest mass hybrid star 
has the radius 10.86 km, the mass 1.81M Q and the cen- 
tral baryon density 7.58po- This star has a large (8 km) 
quark core, and a relatively thin outer layers of hadronic 
matter and a crust. It also appears that the quark cores 
are quite large even for stars with relatively low central 
energy densities. Clearly, this is the reflection of the fact 
that the mixed quark phase starts to develop at rather 



low densities, ps ~ 2.75po7 m the NJL model used. 

The quark core in a hybrid star is separated from 
the lower density hadron-quark layer by a sharp inter- 
face where the baryon density and the energy density 
change by a jump. In a real star, this interface should 
be naturally smoothed over the distances comparable to 
the physical sizes of different matter components in the 
mixed phase. Our estimates show that a typical thick- 
ness of the layers of different quark phases is about 10 
fm. This should also be comparable to the physical sizes 
of components in the hadron-quark mixed phase around 
the interface [2^, |2^. This value sets the scale of the 
actual size of the interface. A more detailed study of the 
interface between the mixed phases is outside the scope 
of this paper. 

In our analysis, we used a rather simple quark model 
and a specific hadronic model. This may suggest that 
most of our conclusions in this paper are not quite gen- 
eral. We would like to argue, however, that is not the 
case. In fact, the main results should remain qualita- 
tively correct and model independent, because they are 
based on a simple and clear physical picture. For exam- 
ple, the fact that the most favorable 2SC quark phase 
appears to be positively charged is related to the nature 
of the phase itself rather than to any specific properties 
of the model. It is this fact that was used previously to 
justify that the CFL phase is more favorable than the 
homogeneous 2SC (or 2SC+s) phase in compact stars. 
Here we turn the same arguments around, and find that 
the mixed phase of normal and 2SC quark matter is the 
most favorable nonstrange globally neutral quark matter. 
Indeed, when converting all quark matter into a color su- 
perconductor is impossible, the next best possibility is to 
convert at least a fraction of it. 

It is fare to mention that, in this paper, we made a 
rather conservative assumption that the strange quark 
is very heavy and it does not appear in quark matter 
even at rather high densities. In particular, the highest 
density reached in our model of the hybrid star without 
causing a gravitational instability is 7.58po. We appre- 
ciate that our assumption may not be justified at such 
large densities. In fact, model calculations of Ref. 0] 
suggest that the strange quarks start to appear in quark 
matter at baryon densities around bpo (this corresponds 
to p = pb/3 ~ 450 MeV). If we accept this, the hybrid 
stars which are heavier than "O-stars" in Figs. [S]and[7| 
should contain some strange matter at their most central 
regions. Apparently, the appearance of the strangeness 
would result in reducing the maximum value of the mass 
of the hybrid star. As one can see from Fig. \7\ however, 
there is no much room for a big change of the maximum 
mass. Of course, adding strangeness may also open the 
possibility of a n ew g eneration of stars with considerably 
smaller radii piEtf! 

In the future, it would be interesting to study the phys- 
ical properties of the hybrid stars that have not been 
addressed in this paper. For example, it would be inter- 
esting to study the neutrino emissivity and the mean free 



11 



path of neutrinos inside quark cores of such stars. The 
corresponding results would be crucial for understanding 
the cooling mechanism of the hybrid stars. Also, it would 
be interesting to address the magnetic and various trans- 
port properties. The nature of the quark mixed phase 
suggests that a flux of the magnetic field should pene- 
trate through the quark core. However, it might also be 
interesting to see how the flux is distributed inside the 
mixed phase, and how this may affect the star properties. 
Finally, it would be useful to generalize the current study 
to the case of rotating stars. 
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